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Abstract. We introduce an approximation scheme to perform an analytic study of the 
oscillation phenomena in a pedagogical and comprehensive way. By using Gaussian wave 
packets, we show that the oscillation is bounded by a time-dependent vanishing func- 
tion which characterizes the slippage between the mass-eigenstate wave packets. We also 
demonstrate that the wave packet spreading represents a secondary effect which play a 
significant role only in the non-relativistic limit. In our analysis, we note the presence of 
a new time-dependent phase and calculate how this additional term modifies the oscillat- 
ing character of the flavor conversion formula. Finally, by considering Box and Sine wave 
packets we study how the choice of different functions to describe the particle localization 
changes the oscillation probability. 

PACS. 02.30. Mv - 03.65.-W 

I. INTRODUCTION 



Recently the great interest in the quantum oscillation phenomena [1-3] has stirred up an increasing 
number of works devoted to several theoretical approaches to particle mixing and oscillations [4~ 
6]. Notwithstanding the exceptional ferment in this field, the conceptual difficulties hidden in the 
oscillation formulas represent an intriguing, and sometimes embarrassing, challenge for physicists. 

The standard plane wave treatment [7,8] is the most elementary approach used to study the flavor 
oscillation problem. However, despite being physically intuitive and simple, it is, strictly speaking, 
neither rigorous nor sufficient for a complete understanding of the physics involved in quantum oscil- 
lations. The plane wave approach implies a perfectly well-known energy-momentum and an infinite 
uncertainty on the space-time localization of the oscillating particle. Oscillations are destroyed under 
these assumptions [9]. In order to overcome such difficulties, an intermediate wave packet model for 
ultra-relavistic neutrinos was introduced by Kayser [9] and followed by other authors [2,6, 10]. Mean- 
while, a common argument against this approach is that oscillating particles are not, and cannot be, 
directly observed [11]. It would be more convincing to write a transition probability between observable 
particles involved in the production and detection of the oscillating particle in an external wave packet 
framework [3,12]. The particle to be studied is represented by a relativistic propagator; it propagates 
between a source and a detector, where wave packets representing the external particles are in inter- 
action. The function which represents the overlap of the incoming and outgoing wave packets in the 
external wave packet model corresponds to the wave function of the propagating mass-eigenstate in 
the intermediate wave packet formalism. Reixiarkably, it could be shown that the probability densities 
for ultra-relativistic stable oscillating particles in both frameworks are mathematically equivalent [3] . 
Thus, it makes sense, in the external wave packet framework, to consider a wave packet associated 
with the propagating particle. However, this wave packet picture brings up a problem, as the overlap 
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function takes into account not only the properties of the source, but also of the detector. This is 
unusual for a wave packet interpretation and not satisfying for causality [3] . This point was clarified by 
Giunti [16] who solves this problem by proposing an improved version of the intermediate wave packet 
model where the wave packet of the oscillating particle is explicitly computed with field-theoretical 
methods in terms of external wave packets. Despite not being applied in a completely free way, the 
(intermediate) wave packet treatment commonly simplifies the discussion of some physical aspects 
going with the oscillation phenomena [13-15]. In this context, we just establish a condensed scheme 
to analytically study the flavor oscillation phenomena, since, in the literature, numerous prescriptions 
are somewhat confusing. 

Quite generally, the analytical approaches for the mass-eigenstate time evolution do not concern 
with the wave packet limitations. In particular, Gaussian wave packets [6, 16] enable us to quantify the 
first and the second order corrections to the oscillation character of propagating particles. In section 
II, we introduce Gaussian wave packets and assume a sharply peaked momentum distribution. Then, 
we approximate the mass-eigenstate energy in order to analytically obtain the expressions for the wave 
packet time evolution and for the flavor oscillation probability. The energy expansion is taken up to 
the second order term and the wave packet spreading and slippage effects are quantified in both non- 
relativistic (NR) and ultra-relativistic (UR) propagation regimes. We also identify an additional time- 
dependent phase which changes the standard oscillating character of the flavor conversion formula. In 
section III, we introduce Box and smoothly vanishing Sine wave packets and study how the choice of 
different function in describing the particle localization could play a significant role in the oscillation 
probability. We draw our conclusions in section IV. 



The main aspects of oscillation phenomena can be understood by studying the two flavor problem. In 
addition, substantial mathematical simplifications result from the assumption that the space depen- 
dence of wave functions is one-dimensional (z-axis). Therefore, we shall use these simplifications to 
calculate the oscillation probabilities. In this context, the time evolution of fiavor wave packets can 
be described by 



where Ua and Ufj are flavor-cigcnstatcs and Ui and 1^2 arc mass-cigenstatcs. The probability of finding 
a flavor state Ufj at the instant t is equal to the integrated squared modulus of the Ufj coefficient 



II. GAUSSIAN WAVE PACKETS 



$(z, t) = (j)i{z, t) cosOvi + (j)2{z, t) sin 6* 1^2 

= \4'i{z,t)cos^ 9 + (j)2{z,t)siv? 9\ Ua + [(j)i{z,t) — (j)2{z,i)]cos6svti9vii 
= (j)a(,z,t;6)ua + (j)p{z,t;6)up, 



(1) 





where iNT(t) represents the mass-eigenstate interference term given by 




(3) 



Let us consider mass-eigenstate wave packets given at time t = by 




s = 1, 2. 



(4) 



The wave functions which describe their time evolution are 





where 




and 
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In order to obtain the oscillation probability, we can calculate the interference term iNT(t) by solving 
the following integral 



TT^ 'fiiPz - Pi)'fi{Pz - P2) exp [-i AE{p^) t] = 



exp 



(aAp)2 



TT^ v'^ [Pz - Po) exp [-i AE{p^) t] 



(6) 



where we have changed the z-integration into a p^-intcgration and introduced the quantities 



Ap = pi-P2, Po^ ^{pi 



-P2) and AE{p^) = E{p;,,mi) - E{pz,m2) . 



The oscillation term is bounded by the exponential function of a Ap at any instant of time. Under 
this condition wc could never observe a pure flavor-cigenstate. Besides, oscillations are considerably 
suppressed if aAp > 1. A necessary condition to observe oscillations is that aAp <IC 1. This constraint 
can also be expressed by Sp 3> Ap where Sp is the momentum uncertainty of the particle. The overlap 
between the momentum distributions is indeed relevant only for Sp 3> Ap. Consequently, without loss 
of generality, we can assume 



iNT(t) = Re 



dpz 
2ti 



V^{Pz - Po) exp [-i AE{p^) t] 



(7) 



In literature, this equation is often obtained by assuming two mass-eigenstate wave packets described 
by the "same" momentum distribution centered around the average momentum Po- This simplifying 
hypothesis also guarantees the instantaneous creation of a pure flavor-eigenstate i/q, at t = [14], 
hence, in what follows, we shall use this simplification. 



II.A. THE ANALYTICAL APPROACH 

In order to obtain an analytic expression for 0s (z, t) by solving the integral in Eq.©, wc firstly rewrite 
the energy E{pz,ms) as 



E{p^,ms) = Es 



9 
Po 



E^ 



Es[l 



(8) 



where 



E,^iml+pl)^, Vs 



Po_ 

E, 



and (Ts 



Pz - Po 
E. 



By assuming a sharply peaked momentum distribution, i. e. (ai^s)^^ ~ ct^ ^ 1, we can expand the 
energy E(pz,ms) in a power series of as- Meanwhile, the integral in Eq.ljSJ can be analytically solved 
only if we consider terms up to order ct^ in the series expansion. In this case, the energy E(pz,ms) is 
approximated by 



E{Pz 



E, 



Es 



Po<Js 



2Es 



(9) 



The zero-order term in the previous expansion, Es, gives the standard plane wave oscillation phase. 
The first-order term, PoCfs , will be responsible for the slippage due to the different group velocities of the 

mass-eigenstate wave packets and represents a linear correction to the standard oscillation phase [14]. 

2 

Finally, the second-order term, ^^c^, which is a (quadratic) secondary correction will give the well- 
known spreading effects in the time propagation of the wave packet and will be also responsible for a 
"new" additional phase to be computed in the final calculation. In the case of Gaussian momentum 
distributions for the mass-eigenstate wave packets, these terms can all be analytically quantified. 
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By substituting © in Eq.Q and changing the p^-integration into a (Ts-integration, we obtain the 
expHcit form of the mass-eigenstate wave packet time evolution, 



{z, t) K. {2tt a^) 4 exp [-i{Es t - po z)] 



where 



Eg exp 



,2 tp2 „2 
^ Es <^s 



exp 



-i {po t - Es z)<7s - i TT-^cr^ 



2Eg 



exp [— z [Eg t — Po z)] exp 



[z - V, tf 



iOs(t,z) 



(10) 



4 m 

as{t) ^ a ^ 1 -\ — ] ^i^d Og{t, z) = ^ - arctan 



2mit 



a2 £3 



a' El am 



The timc-dcpcndcnt quantities as{t) and 6s{t,z) contain all the physically significant information 
which arise from the second order term in the power series expansion The spreading of the 
propagating wave packet can be immediately quantified by interpreting as(t) as a time-dependent 
width, i. e. the spatial localization of the propagating particle is effectively given by as{t) which 
increases during the time evolution. In the NR propagation regime, as{t) is approximated by a""(t) = 

ay^l + ^^4^^^ [19]- For times t ^ a^mg the effective wave packet width a^'^{t) becomes much larger 

than the initial width a. Otherwise, the wave packet spreading in the UR propagation regime is 

is approximated by a""'{t) = ay^lH^"^^^ w a. The UR spreading is practically negligible if we 

consider the same time-scale T for both NR and UR cases, i. e. a""-{T) ^ a""(T). To illustrate this 
characteristic, we plot the time-dependence of ag{t) in Fig^ where we have assumed a particle with 
a definite mass value irig. By computing the squared modulus of the mass-eigenstate wave function. 



exp 



2{z~Yg 



(11) 



we illustrate the wave packet spreading in both NR and UR propagation regimes in Fig|3 which is in 
correspondence with Fig^ It confirms that the wave packet spreading is irrelevant for UR particles. 

Returning to Eq. (|10|) , we could interpret another second order effect by observing the time-behavior 
of the phase Os(t, z). By taking into account the wave packet localization, we assume that the amplitude 
of the wave function is relevant in the interval \z — Vg i| < as{t). Due to the z-dependence, each wave 
packet space-point z evolves in time in a different way. If we observe the propagation of the space- 
point z = Vs i, the crescent function 6g{t, Vgt) assume values limited by the interval [0, j[. Otherwise, 
for any other space-point given by z = v^ t -I- K as{t), < \K\ < 1, the phase 9s{t, z) does not have 
a lower limit. We shall show in the next subsection that the presence of a time-dependent phase can 
modify the oscillation character of the flavor conversion formula. Anyway, the phase 0g{t,z) is not 
influent on the free mass-eigenstate wave packet propagation as we can see from Ea. (|ll|l . 



II. B. THE OSCILLATION PROBABILITY 

After having analytically quantified the second order corrections to the time evolving mass-eigenstate 
wave packets, we now compute the interference term iNT(t) in order to obtain an explicit expression 
for the flavor conversion probability. By solving the integral with the approximation (jHJ and 
performing some mathematical manipulations, wc obtain 

iNT(t) = BND(t) X OSC(t), (12) 

where we have factored the time- vanishing bound of the interference term given by 

BND(i) = [1 + Sp^(t)] exp 



(Avt)2 



2a2 [1 + Sp2(t)] 



(13) 
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and the time-oscillating character of the flavor conversion formula given by 

Osc(<) =i?e{exp [-iAEt-iQ{t)]} = cos [AEt + e{t)] 

where 



r. / X t ^ f Avi 



and Q{t) = 



Sp3(0 



1 a^v^ 
- arctan [Sp(i)] - -t-t^t 5 — t 

2 ^ ^ 2p2 1 + Sp2 t) 



with 



AE\ 



i- and 



£■ — \/ EiE2- 



(14) 



(15) 



(16) 



The time-dependent quantities Sp(t) and 9(t) carry the second order corrections and, consequently, 
the spreading effect to the oscillation probability formula. If Ai? ^ E, the parameter p is limited by 

the interval [1,-2] and it assumes the zero value when w i. Therefore, by considering increasing 

values of po, from NR to UR propagation regimes, and fixing the time derivatives of Sp(t) and 

Q{t) have their signals inverted when ^ reaches the value ^. 

To simplify our presentation, let us study separately the time-dependent functions BND(t) and 
Osc(t). The slippage between the mass-eigenstate wave packets is quantified by the vanishing behavior 
of BND(t). 

In order to compare BND(t) with the correspondent function without the second order corrections 
(without spreading), 

\ (Avi)2" 

BND„,s(t) = CXp 



2a2 

we substitute Sp(i) given by the expression H14|l in Ea. ljl^l) and we obtain the ratio 



BND(i) 

BnDv^s(<) 



AEt 
a2 E^ 



exp 



2a6 E» 



1 



(17) 



(18) 



The NR limit is obtained by setting = I and po = in Eq. 118|) . In the same way, the UR limit is 
obtained by setting = 4 and Po = E. 

In fact, the minimal influence due to second order corrections occurs when ^ « i (p « 0). 
Returning to the exponential term of Ea. ljlHII . we observe that the oscillation amplitude is more 
relevant when Avt <^ a. It characterizes the minimal slippage between the mass-eigenstate wave 
packets which occur when the complete spatial intersection between themselves starts to diminish 



BND(t) 



ws{t) 



1. 



during the time evolution. Anyway, under minimal slippage conditions, we always have 

We plot the ratio given in Eq. (|18|l for different propagation regimes in Fig|31 where we have ar- 
bitrarily set aE = 10. For asymptotic times, the time-dependent term Sp(t) effectively extends the 
interference between the mass-eigenstate wave packets since 



aE 



BND(i) t-*a 

BND„,s(t) ^ {pAEty. 



■ exp 



2a2 £4 



» 1, 



(19) 



but, in this case, the oscillations are almost completely destroyed by BND(i) (see Fig.©). 

The oscillating function Osc(i) of the interference term iNT(t) differs from the standard oscillating 
term, cos[AEt], by the presence of the additional phase Q{t) which is essentially a second order 
correction. The modifications introduced by the additional phase Q{t) are presented in Fig01 where 
we have compared the time-behavior of Osc(t) to cos[A£'t] for different propagation regimes. To 
study the phase Q{t), let us conveniently define a time t = to > which sets the zero of &{t), i. e. 
&{to) = 0. If t < to, the modulus of the phase 8(t) reaches an upper limit when 
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therefore, the maxhnum of |6(i)| depends, not only on the propagation regime {po value), but also on 
the wave packet width a. Anyway, the values assumed by \Q{t)\ are restricted to the interval [0, f [. 
Otherwise, if t > to, the phase 0(<) does not have a limit and its time-dependence is essentially 
given by the second term of Ea. (|15|l . However, it is important to notice that for t > to the oscillating 
character is gradually destroyed by BND(t). Consequently, another bound effective value assumed by 
&{t) is determined by the vanishing behavior of BND(t). To illustrate this point, we plot both the 
curves representing BND(t) and Q{t) in Fig[Slby considering the same parameters used in the study of 
BND(t). We note the phase slowly changing in the NR regime. The modulus of the phase |6(i)| rapidly 

reaches its upper limit when ^ > 5 and, after a time t = to, it continues to evolve approximately 
linearly in time. But, effectively, the oscillations rapidly vanishes after t = to- 

By superposing the effects of BND(t) and the oscillating character Osc(t) expressed in Fig[Sl we 
immediately obtain the flavor oscillation probability which is explicitly given by 



sin^ [26] 



1 - [1 + Sp^{t)] * exp 



(Avi)2 



2a2 [1 + Sp2(i)] 



cos [AEt + Q{t)] } .(21) 



Obviously, the larger is the value of aS, the smaller are the wave packet effects. If it was sufR- 
ciently larger to not consider the second order corrections expressed in Eq.®, we could compute the 
oscillation probability with the leading corrections due to the slippage effect. 



sin^ \2e] 



1 — exp 



(Avt)2 



2a2 



cos [ Ai; t] 



(22) 



which corresponds to the same result obtained by [14]. Under minimal slippage conditions, the above 
expression reproduces the standard plane wave result. 



, sin^ \2e] , 



{Avtf 
2a2 



sin^ [2611 r r, n 

—LA{l-cos[AEt]}, 



cos \AE t] 



(23) 



since we have assumed a E ^ 1. 



III. ANALYSIS WITH DIFFERENT WAVE PACKETS 

In this section we verify in what circumstances the form of the wave function can change the flavor 
oscillation probability. To describe the wave packet time evolution, let us now consider a Box function 
and a (smoothly vanishing) Sine function in the place of a Gaussian function. In the previous section, 
we have noticed it is remarkably simple to perform an analytical study with a Gaussian wave packet 
since its Fourier transform (FT) in the momentum space is also a Gaussian function. In opposition, the 
analytical study with Box and Sine functions constrain us to perform the calculations by considering 
only the first order corrections in Eq. ©, i. c. 

E{p,^,ms)Ki Es+pods (24) 

which only sets the slippage leading term. We can observe from Fig|3 that considering only the 

first order corrections results in a good approximation for propagation regimes where -g^ > 5 since 
the oscillations are almost completely destroyed after any relevant second order correction takes 

place. Besides, for NR propagation regimes, i e. when ^ < |, by observing the Figl^J we have 
already noticed that the first and second order approximations are equivalent under minimal slippage 

conditions f-n^^^^^^ « iV 

To simplify the discussion, we shall adopt the following definition for the initial state. 



0«(z,O) = ^^«(z) exp[zp„z], 



(25) 
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where i = G, B, S correspond respectively to Gaussian, Box and Sine functions. The wave packet time 

evolution will be expressed in terms of (p^'''^{pz ~ Po) which is the FT of 4>^s \z, 0), and the oscillation 
probability will be immediately computed through the expression 

As we have seen in the previous section, in the case of a Gaussian function, we have 



f 2 V 




2 









— 

and ^^^''ipz - Po) = (2^0) " 6xp 



al {pz - Pof 



In this case, the wave packet has the form 

1 

2 



11 at 



exp [— z [Es t ~ Po z)] exp 



[Z - Vs tf 



(26) 



and the oscillation probability is reproduced by Ea. l|22l) . Obviously, such results could be directly 
obtained by setting as{t) = a and Os{t, z) = in Ea. Hl()|l . 
In the case of a Box function wc have 



z e 



flB as 
2 ' 2 . 



and ip^'''>{pz - Po) 







ai ipz - Po) 



Ab {Pz -Po) 



In this case, the wave packet has the form 

Ob ^ exp [-i {Es t~po z)] 



4''Hz,t) 



and the oscillation probability becomes 



G [vst- 



(27) 







Ob. 
2 . 



sin^ [29 



1 - 



Finally, in the case of a Sine function we have 



I — Avt 

Ob 



cos[A£'t]| i<|| 



(28) 



^ ^ AV 



as \ 2 sm z Cc, 



and .^^^^(p^-po) 




as iPz-Po) e [-1, 1] 
fls iPz-Po) i [-1, 1] 



In this case, the wave packet has the form 



ct>'^^\z,t) 



exp [-i {Es t- poz] 



sin [og (z — Vs t)\ 



(z - Vs i) 



and the oscillation probability becomes 

sin^ \2e] 



P^^\u^^vp-t) 



1 - 



Avi 



sm 



Avi 

as 



cos [Ai; t] 



(29) 



(30) 



The above results deserve some comments. Firstly, we observe that all the three wave packet forms 
give the same oscillating character. In a simplified analysis, independently of the propagation regime 
and without setting any parameter value, we can compare the vanishing character of each oscillation 
probability in terms of a common variable x{t) = Ail^ gy (defining the coefficients q;b = ^ and 



as 



and recovering the definition of Bnd(<) we can write 



BND''^^(t) = exp 



x^{t) 



, BND(^'(t) 



1 — a^x{t) a^x{t) < 1 







,{t) > 1 



and BND''''(i) 



sm \ a<,x 



it)] 



as x{t) 



8 



Alex E. Bernardini, Stefano De Leo: An analytic approach to the wave packet formalism 



Under minimal slippage conditions, i. c. when x{t) ^ 1, Bnd''°-'(<:) and BND'-^-'(t) vanish quadratically. 
Particularly, if we had set = a/S, we would have 

Bnd^'^'CO = BND("'(i) « 1 - (31) 

i. e. under minimal slippage conditions, Gaussian and Sine functions would give exactly the same 
oscillation probabilities. To summarize the above results, we show the oscillation probabilities by 
considering the three wave packet forms in Fig|Sl where we have adopted = 1 and as = V3- 
Predominantly for Sine functions, there will always be a reminiscent oscillating character during the 
particle propagation. In opposition, BND(t)'^'(t) vanishes linearly and the correspondent oscillation 
probability goes much more rapidly to zero. Its oscillating character is suddenly ended up when x{t) = 
The Sine wave packets still provide another peculiar behavior. Their correspondent oscillations 
vanish at each zero of sin [x(t)] but the probability returns to oscillate. After each intermediary zero, 
the function sin [x(t)] changes the signal itself, consequently, its maximum and minimum values are 
interchanged. In Fig[7| we illustrate the correspondent slippage between the mass-eigenstate wave 
packets for each case. 



IV. CONCLUSIONS 

In this paper we have analytically computed the second order modifications to the flavor conversion 
formula by using Gaussian wave packets. Under the particular assumption of a sharply peaked momen- 
tum distribution, we have obtained an explicit expression for the time evolution of the mass-eigenstates 
and identified the wave packet spreading for (U)R and NR propagation regimes. In particular, we have 
observed that the spreading represents a minor modification effect which is practically irrelevant for 
(ultra)relativistic propagating particles. We have also observed the presence of an additional time- 
dependent phase in the oscillating term of the flavor conversion formula. Such an additional phase 
presents an analytic dependence on time which changes the oscillating character in a peculiar way. 
These modifications are less relevant when p^ sa ^E'^ and more relevant for NR propagation regimes 
Anyway, they become completely irrelevant for UR propagatiton regimes due to the vanishing behav- 
ior of the interference term in the oscillation probability formula. Some influences of this additional 
phase on the oscillation problem were already appointed in reference [17]. 

We know, however, that our results are strongly influenced by the Gaussian wave packet choice. In 
order to understand how the wave packet form modifies the oscillation probability, we have quantified 
the slippage between the mass-eigenstate wave packets by studying a Box and a Sine localization. 
In fact, by following a first order analytic approximation, a simple comparison among the different 
vanishing character of the oscillation probability formulas has illustrated that, under minimal slippage 
conditions, the Sine and the Gaussian functions provide similar results whereas the Box function makes 
the oscillations vanishing more rapidly. 

To conclude, we emphasize that, an analytical study complements and clears up several aspects 
already introduced in the study of quantum oscillation phenomena. 
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FIGURES 
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Fig. 1. The time-dependence of the wave packet width as{t) is given for different values of the ratio po / nis. 
By considering a fixed mass value nis, we compare the non-relativistic {po <^ rUs) and tire ultra-relativistic 
{Po 3> rUs) propagation regimes. We observe that the spreading is much more relevant in the former case. In 
the ultra-relativistic limit {nis = 0), the wave packet does not spread and as{t) assumes a constant value a. 
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Fig. 2. The wave packet spreading in both non-relativistic and ultra-relativistic propagation regimes is de- 
scribed at time t = Aa^rris in correspondence with FigQ The solid line represent the shape of the wave packet 

2 

at time t — 0. In the case of an ultra-relativistic propagation expressed in terms of = 10, the spreaa 
indeed irrelevant. 
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Fig. 3. The comparison between the vanishing behavior with (BND(t)) and without (BNDivs(t)) the second 
order corrections for different propagation regimes. In order to have a realistic interpretation of the information 
carried by the second order corrections we arbitrarily fix a_E = 10. The second order corrections could be 
indeed effective for both non-relativistic and (ultra)relativistic propagation regimes, however, the oscillations 

are destroyed much more rapidly in the latter case. 



Tf 



i, the second order corrections are minimal. 
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Fig. 4. The time-behavior of Osc(t) compared with the standard plane wave oscillation given by cos [AEt] 
for different propagation regimes. The additional phase &{t) changes the oscillating character after some time 

of propagation. The maximal deviation occurs for ^ ~ § . 
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Fig. 5. The time-behavior of the additional phase 0(t). The values assumed by Q{t) are effective while the 
interference term does not vanish. In the upper box we can observe the behavior of BND(t) which determines the 

2 

limit values effectively assumed by 0(t) for each propagation regime. For relativistic regimes with > 5, the 
function 0(t) rapidly reaches its lower limit as we can observe in the small box above. We have used aE = 10. 
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Fig. 6. The flavor conversion probabilities for Gaussian, Box and Sine wave packets by taking into account the 
first order correction in an analytical calculation of iNT(t). By assuming ac = as = -^as, the Gaussian and 
the Sine wave packets provide exactly the same quadratic time dependence under minimal slippage conditions 
whereas the Box wave packets give a completely different behavior where the oscillation probability vanishes 
much more rapidly. 
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Slippage between mass-eigenstate wave packets 



Fig. 7. The slippage between Gaussian, Box and Sme wave packets. We can observe that the interference 
between the Box wave packets is abruptly interrupted while the other two wave packets continue to interfere 
during longer times. It completes the explanation of the oscillation behavior illustrated in Fig|S] 



